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^ . Abstract 

The highest energy cosmic rays observed possess macroscopic energies and their origin is hkely to 
be associated with the most energetic processes in the Universe. Their existence triggered a flurry of 
1/-^ . theoretical explanations ranging from conventional shock acceleration to particle physics beyond the 
I Standard Model and processes taking place at the earliest moments of our Universe. Furthermore, 
many new experimental activities promise a strong increase of statistics at the highest energies and 
. a combination with 7— ray and neutrino astrophysics will put strong constraints on these theoretical 
I models. We give an overview over this quickly evolving research field with focus on testing new particle 
physics. 

p. 



^ '. 1 Introduction 



Over the last few years, several giant air showers have been detected confirming the arrival of cosmic 
^ . rays (CRs) with energies up to a few hundred EeV (1 EeV = 10^® eV) |2|, ^. The existence of 
^ I such ultra-high energy cosmic rays (UHECRs) pose a serious challenge for conventional theories of CR 
origin based on acceleration of charged particles in powerful astrophysical objects. The question of 
the origin of these UHECRs is, therefore, currently a subject of much intense debate and discussions 
as well as experimental efforts; see Ref. [Q] for recent brief reviews, and Ref. |||] for a detailed review. 

The problems encountered in trying to explain UHECRs in terms of "bottom-up" acceleration 
mechanisms have been well-documented in a number of studies; see, e.g., Refs. j^, |^, |8|. It is hard 
to accelerate protons and heavy nuclei up to such energies even in the most powerful astrophysical 
objects such as radio galaxies and active galactic nuclei. Also, nucleons above ~ 70 EeV lose energy 
drastically due to photo-pion production on the cosmic microwave background (CMB) — the Greisen- 
Zatsepin-Kuzmin (GZK) effect Q — which limits the distance to possible sources to less than ~ 
100 Mpc [^. Heavy nuclei are photodisintegrated in the CMB within a few Mpc |1C]. There are no 
obvious astronomical sources within ~ 100 Mpc of the Earth |11, |^]. 



The distance restriction imposed by the GZK effect can be circumvented if the problem of energetics 
is somehow solved separately and if one postulates new particles beyond the Standard Model; this 
will be discussed in the third section. 

In contrast, in the "top-down" scenarios, which will be discussed in the last section, the problem of 
energetics is trivially solved. Here, the UHECR particles are the decay products of some supermassive 
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"X" particles of mass mx ^ 10^^ eV, and have energies all the way up to ~ nix- Thus, no acceleration 
mechanism is needed. The massive X particles could be metastable relics of the early Universe with 
lifetimes of the order the current age of the Universe or could be released from topological defects that 
were produced in the early Universe during symmetry-breaking phase transitions envisaged in Grand 
Unified Theories (GUTs). If the X particles themselves or their sources cluster similar to dark matter, 
the dominant observable UHECR contribution would come from the Galactic Halo and absorption 
would be negligible. 

The main problem of non-astrophysical solutions of the UHECR problem in general is that they 
are highly model dependent. On the other hand, they allow to at least test new physics beyond the 
Standard Model of particle physics (such as Grand Unification and new interactions beyond the reach 
of terrestrial accelerators) as well as early Universe cosmology (such as topological defects and/or 
massive particle production in inflation) at energies often inaccessible to accelerator experiments. 

The physics and astrophysics of UHECRs are intimately linked with the emerging field of neutrino 



astronomy (for reviews see Refs. [12, |T^) as well as with the already established field of 7— ray 
astronomy (for reviews see, e.g., Ref. [0). Indeed, all scenarios of UHECR origin, including the 
top-down models, are severely constrained by neutrino and 7— ray observations and limits. In turn, 
this linkage has important consequences for theoretical predictions of fluxes of extragalactic neutrinos 
above a TeV or so whose detection is a major goal of next-generation neutrino telescopes: If these 
neutrinos are produced as secondaries of protons accelerated in astrophysical sources and if these 
protons are not absorbed in the sources, but rather contribute to the UHECR flux observed, then 
the energy content in the neutrino flux can not be higher than the one in UHECRs, leading to the 
so called Waxman Bahcall bound for sources with soft acceleration spectra [15, 17 [. If one of these 
assumptions does not apply, such as for acceleration sources with injection spectra harder than E~'^ 
and/or opaque to nucleons, or in the top-down scenarios where X particle decays produce much fewer 
nucleons than 7— rays and neutrinos, the Waxman Bahcall bound does not apply, but the neutrino 
flux is still constrained by the observed diffuse 7— ray flux in the GeV range. 



2 Propagation of Ultra-High Energy Radiation: Simu- 
lations 

Before discussing specific scenarios for UHECR origin we give a short account of the numerical tools 



used to compute spectra of ultra-high energy cosmic and 7— rays, and neutrinos [p^ , 19, |M. In 
the following we assume a flat Universe with a Hubble constant of -ff = 70 km sec~^Mpc~ and a 
cosmological constant = 0.7, as favored by current observations. 

The relevant nucleon interactions implemented are pair production by protons {p^h — > pe~e^), 
photoproduction of single or multiple pions (A^7fe N nvr, n > 1), and neutron decay. 

7-rays and electrons/positrons initiate electromagnetic (EM) cascades on low energy radiation 
fields such as the CMB. The high energy photons undergo electron-positron pair production (PP; 
776 e~e^), and at energies below ~ 10^^ eV they interact mainly with the universal infrared and 
optical (IR/0) backgrounds, while above ~ 100 EeV they interact mainly with the universal radio 
background (URB). In the Klein-Nishina regime, where the CM energy is large compared to the 
electron mass, one of the outgoing particles usually carries most of the initial energy. This "leading" 
electron (positron) in turn can transfer almost all of its energy to a background photon via inverse 
Compton scattering (ICS; 67;, — > e'7). EM cascades are driven by this cycle of PP and ICS. The energy 
degradation of the "leading" particle in this cycle is slow, whereas the total number of particles grows 
exponentially with time. All EM interactions that influence the 7-ray spectrum in the energy range 
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10^ eV < E < 10^^ eV, namely PP, ICS, triplet pair production (TPP; e^b ee~e^), and double 
pair production (DPP, 77;, — > e~e^e~e'^), as well as synchrotron losses of electrons in the large scale 
extragalactic magnetic field (EGMF), are included. 

Similarly to photons, UHE neutrinos give rise to neutrino cascades in the primordial neutrino 
background via exchange of W and Z bosons |21]. Besides the secondary neutrinos which drive the 
neutrino cascade, the W and Z decay products include charged leptons and quarks which in turn 
feed into the EM and hadronic channels. Neutrino interactions become especially significant if the 
relic neutrinos have masses rriu in the eV range and thus constitute hot dark matter, because the Z 
boson resonance then occurs at an UHE neutrino energy Sres = 4 x 10^^(eV/m,y) eV. In fact, the 
decay products of this "Z-burst" have been proposed as a significant source of UHECRs ||2^. The big 
drawback of this scenario is the need of enormous primary neutrino fluxes that cannot be produced by 
known astrophysical acceleration sources [|l9[, and thus most likely requires a more exotic top-down 
type source such as X particles exclusively decaying into neutrinos [^^. Even this possibility appears 
close to being ruled out due to a tendency to overproduce the diffuse GeV 7— ray flux observed by 
EGRET H, H 

The two major uncertainties in the particle transport are the intensity and spectrum of the URB 



for which there exist no direct measurements in the relevant MHz regime [25, |23], and the average 



value of the EGMF. Simulations have been performed for different assumptions on these. A strong 
URB tends to suppress the UHE 7-ray flux by direct absorption whereas a strong EGMF blocks EM 
cascading (which otherwise develops efficiently especially in a low URB) by synchrotron cooling of the 
electrons. For the IR/0 background we used the most recent data [p^. 

In top-down scenarios, the particle injection spectrum is generally dominated by the "primary" 
7-rays and neutrinos over nucleons. These primary 7-rays and neutrinos are produced by the decay 
of the primary pions resulting from the hadronization of quarks that come from the decay of the X 
particles. In contrast, in acceleration scenarios the primaries are accelerated protons or nuclei, and 
7-rays, electrons, and neutrinos are produced as secondaries from decaying pions that are in turn 
produced by the interactions of nucleons with the CMB. 



3 New Primary Particles and Interactions 

A possible way around the problem of missing counterparts within acceleration scenarios is to propose 
primary particles whose range is not limited by interactions with the CMB. Within the Standard 
Model the only candidate is the neutrino, whereas in extensions of the Standard Model one could 
think of new neutrals such as axions or stable supersymmetric elementary particles. Such options are 
mostly ruled out by the tension between enforcing small EM coupling and large hadronic coupling to 
ensure normal air showers [^^. Also suggested have been new neutral hadronic bound states of light 
gluinos with quarks and gluons, so-called R-hadrons that are heavier than nucleons, and therefore 



have a higher GZK threshold |2^. Since this too seems to be disfavored by accelerator constraints [30| 
we will here focus on neutrinos. 

In both the neutrino and new neutral stable particle scenario the particle propagating over extra- 
galactic distances would have to be produced as a secondary in interactions of a primary proton that 
is accelerated in a powerful AGN which can, in contrast to the case of extensive air showers (EAS) 
induced by nucleons, nuclei, or 7— rays, be located at high redshift. Consequently, these scenarios 
predict a correlation between primary arrival directions and high redshift sources. In fact, possible ev- 
idence for a correlation of UHECR arrival directions with compact radio quasars and BL-Lac objects, 
some of them possibly too far away to be consistent with the GZK effect, was recently reported [pi] ]. 
Only a few more events could conflrm or rule out the correlation hypothesis. Note, however, that 
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these scenarios require the primary proton to be accelerated up to at least 10^^ eV, demanding a very 
powerful astrophysical accelerator. 

3.1 New Neutrino Interactions 

Neutrino primaries have the advantage of being well established particles. However, within the 
Standard Model their interaction cross section with nucleons, whose charged current part can be 
parametrized by 

4n{E) ~ 2.36 X 10-=^2(^/ioi9 eV)°-363 cm^ , (1) 

for 10^^ eV ^ E < 10^^ eV, falls short by about five orders of magnitude to produce ordinary air 
showers. However, it has been suggested that the neutrino-nucleon cross section, a^N, can be enhanced 
by new physics beyond the electroweak scale in the center of mass (CM) frame, or above about a PeV 
in the nucleon rest frame. Neutrino induced air showers may therefore rather directly probe new 
physics beyond the electroweak scale. 

One possibility consists of a large increase in the number of degrees of freedom above the elec- 
troweak scale [^]. A specific implementation of this idea is given in theories with n additional large 
compact dimensions and a quantum gravity scale M^^n ~ TeV that has recently received much at- 
tention in the literature because it provides an alternative solution (i.e., without supersymmetry) 
to the hierarchy problem in grand unifications of gauge interactions. It turns out that the largest 
contribution to the neutrino-nucleon cross section is provided by the production of microscopic black 
holes centered on the brane representing our world, but extending into the extra dimensions. The 
cross sections can be larger than the Standard model one by up to a factor 100 [35|. This is not 
sufficient to explain the observed UHECR events |36]. 



However, the UHECR data can be used to put constraints on cross sections satisfying a^N{E > 
10^^ eV) < 10~'^''cm^. Particles with such cross sections would give rise to horizontal air showers 
which have not yet been observed. Resulting upper limits on their fiuxes assuming the Standard 
Model cross section Eq. (|l|) are shown in Fig. |^. Comparison with the "cosmogenic" neutrino fiux 
produced by UHECRs interacting with the CMB then results in upper limits on the cross section which 
are about a factor 1000 larger than Eq. (|l|) in the energy range between ~ 10^'' eV and ~ 10^^ eV p^]. 
The projected sensitivity of future experiments shown in Fig. § below indicate that these limits could 
be lowered down to the Standard Model one |43|. In case of a detection of penetrating events the 
degeneracy of the cross section with the unknown fiux could be broken by comparing horizontal air 



showers with Earth skimming events |44|. 



4 Top-Down Scenarios 

As mentioned in the introduction, all top-down scenarios involve the decay of X particles of mass 
close to the GUT scale which can basically be produced in two ways: If they are very short lived, as 
usually expected in many GUTs, they have to be produced continuously. The only way this can be 
achieved is by emission from topological defects left over from cosmological phase transitions that may 
have occurred in the early Universe at temperatures close to the GUT scale, possibly during reheating 
after infiation. Topological defects necessarily occur between regions that are causally disconnected, 
such that the orientation of the order parameter associated with the phase transition, can not be 
communicated between these regions and consequently will adopt different values. Examples are 
cosmic strings, magnetic monopoles, and domain walls. The defect density is consequently given by 
the particle horizon in the early Universe. The defects are topologically stable, but time dependent 
motion leads to the emission of particles with a mass comparable to the temperature at which the 
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phase transition took place. The associated phase transition can also occur during reheating after 
inflation. 

Alternatively, instead of being released from topological defects, X particles may have been pro- 
duced directly in the early Universe and, due to some unknown symmetries, have a very long lifetime 
comparable to the age of the Universe. In contrast to Weakly-Interacting Massive Particles (WIMPS) 
below a few hundred TeV which are the usual dark matter candidates motivated by, for example, 
supersymmetry and can be produced by thermal freeze out, such superheavy X particles have to be 
produced non-thermally (see Ref. |^5| for a review). In all these cases, such particles, also called "WIM- 
PZILLAs", would contribute to the dark matter and their decays could still contribute to UHECR 
fluxes today, with an anisotropy pattern that reflects the dark matter distribution in the halo of our 
Galaxy. 

It is interesting to note that one of the prime motivations of the inflationary paradigm was to 
dilute excessive production of "dangerous relics" such as topological defects and superheavy stable 
particles. However, such objects can be produced right after inflation during reheating in cosmologi- 
cally interesting abundances, and with a mass scale roughly given by the inflationary scale which in 
turn is fixed by the CMB anisotropics to ~ 10^^ GeV . The reader will realize that this mass scale 
is somewhat above the highest energies observed in CRs, which implies that the decay products of 
these primordial relics could well have something to do with UHECRs which in turn can probe such 
scenarios! 

For dimensional reasons the spatially averaged X particle injection rate can only depend on the 
mass scale mx and on cosmic time t in the combination 

hx{t) = Km^^t-^+P , (2) 

where k and p are dimensionless constants whose value depend on the specific top-down scenario |^] , 
For example, the case p = 1 is representative of scenarios involving release of X particles from topo- 
logical defects, such as ordinary cosmic strings |^^, necklaces [^] and magnetic monopoles This 
can be easily seen as follows: The energy density ps in a network of defects has to scale roughly 
as the critical density, ps oc pcrit oc where t is cosmic time, otherwise the defects would either 
start to overdose the Universe, or end up having a negligible contribution to the total energy den- 
sity. In order to maintain this scaling, the defect network has to release energy with a rate given 
by Ps = —aps/t oc where a = 1 in the radiation dominated aera, and a = 2/3 during matter 
domination. If most of this energy goes into emission of X particles, then typically k ~ 0{1). In 
the numerical simulations presented below, it was assumed that the X particles are nonrelativistic at 
decay. 

The X particles could be gauge bosons, Higgs bosons, superheavy fermions, etc. depending on 
the specific GUT. They would have a mass mx comparable to the symmetry breaking scale and 
would decay into leptons and/or quarks of roughly comparable energy. The quarks interact strongly 
and hadronize into nucleons (A^s) and pions, the latter decaying in turn into 7-rays, electrons, and 
neutrinos. Given the X particle production rate, dnx/dt, the effective injection spectrum of particle 
species a [a = 7, A, e^,i^) via the hadronic channel can be written as {dnx /dt){2/mx){dNa/dx), 
where x = 2E/mx, and dNa/dx is the relevant fragmentation function (FF). 

We adopt the Local Parton Hadron Duality (LPHD) approximation according to which the 
total hadronic FF, dN^/dx, is taken to be proportional to the spectrum of the partons (quarks/gluons) 
in the parton cascade (which is initiated by the quark through perturbative QCD processes) after 
evolving the parton cascade to a stage where the typical transverse momentum transfer in the QCD 
cascading processes has come down to ~ ^ few hundred MeV, where R is a typical hadron 
size. The parton spectrum is obtained from solutions of the standard QCD evolution equations in 
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modified leading logarithmic approximation (MLLA) which provides good fits to accelerator data at 
LEP energies [^]. Within the LPHD hypothesis, the pions and nucleons after hadronization have 
essentially the same spectrum. The LPHD does not, however, fix the relative abundance of pions and 
nucleons after hadronization. Motivated by accelerator data, we assume the nucleon content /at of the 
hadrons to be ~ 10%, and the rest pions distributed equally among the three charge states. Recent 
work suggests that the nucleon-to-pion ratio may be significantly higher in certain ranges of x values 



at the extremely high energies of interest here |51], but the situation is not completely settled yet. 



4.1 Predicted Fluxes 

Fig. ^ shows results for the time averaged nucleon, 7— ray, and neutrino fluxes in a typical TD scenario, 
along with low energy 7— ray flux constraints and neutrino flux sensitivities of future experiments. The 
spectrum was optimally normalized to allow for an explanation of the observed UHECR events, assum- 
ing their consistency with a nucleon or 7— ray primary. The flux below < 2 x 10^^ eV is presumably due 
to conventional acceleration in astrophysical sources and was not fit. The PP process on the CMB de- 
pletes the photon flux above 100 TeV, and the same process on the IR/0 background causes depletion 
of the photon flux in the range 100 GeV-100 TeV, recycling the absorbed energies to energies below 
100 GeV through EM cascading. The predicted background is not very sensitive to the specific IR/0 
background model, however |52]. The scenario in Fig. || obviously obeys all current constraints within 
the normalization ambiguities and is therefore quite viable. Note that the diffuse 7— ray background 



measured by EGRET |24] up to 10 GeV puts a strong constraint on these scenarios, especially if there 
is already a significant contribution to this background from conventional sources such as unresolved 
7— ray blazars However, this constraint is much weaker for TDs or decaying long lived X particles 
with a non-uniform clustered density ||5^ 

The energy loss and absorption lengths for UHE nucleons and photons are short (< 100 Mpc). 
Thus, their predicted UHE fiuxes are independent of cosmological evolution. The 7— ray flux below 
~ 10^^ eV, however, scales as the total X particle energy release integrated over all redshifts and 
increases with decreasing p []55| . For mx = 2 x 10^^ GeV, scenarios with p < 1 are therefore ruled out, 
whereas constant comoving injection models {p = 2) are well within the limits. 

It is clear from the above discussions that the predicted particle fiuxes in the TD scenario are 
currently uncertain to a large extent due to particle physics uncertainties (e.g., mass and decay modes 
of the X particles, the quark fragmentation function, the nucleon fraction /tv, and so on) as well 
as astrophysical uncertainties (e.g., strengths of the radio and infrared backgrounds, extragalactic 
magnetic fields, etc.). More details on the dependence of the predicted UHE particle spectra and 
composition on these particle physics and astrophysical uncertainties are contained in Ref. |63|. 



We stress here that there are viable TD scenarios which predict nucleon fluxes that are comparable 
to or even higher than the 7— ray fiux at all energies, even though 7— rays dominate at production. 
This occurs, e.g., in the case of high URB and/or for a strong EGMF, and a nucleon fragmentation 
fraction of ~ 10%. Some of these TD scenarios would therefore remain viable even if UHECR induced 
EAS should be proven inconsistent with photon primaries (see, e.g., Ref. [0]). This is in contrast to 
scenarios with decaying massive dark matter in the Galactic halo which, due to the lack of absorption, 
predict compositions directly given by the fragmentation function, i.e. domination by 7— rays. 

The normalization procedure to the UHECR flux described above imposes the constraint Quhecr ~ 
10-22 eV cm ^ sec ^ within a factor of a few ^ ^ for the total energy release rate Qq from TDs 



at the current epoch. In most TD models, because of the unknown values of the parameters involved, 
it is currently not possible to calculate the exact value of Qq from first principles, although it has been 
shown that the required values of Qq (in order to explain the UHECR fiux) mentioned above are quite 
possible for certain kinds of TDs. Some cosmic string simulations and the necklace scenario suggest 
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that defects may lose most of their energy in the form of X particles and estimates of this rate have 
been given 48|. If that is the case, the constraint on Quhecr translates via Eq. (|2|) into a limit on 
the symmetry breaking scale 77 and hence on the mass mx of the X particle: 77 ~ mx ^ 10^'^ GeV []6^ ] . 
Independently of whether or not this scenario explains UHECR, the EGRET measurement of the 
diffuse GeV 7— ray background leads to a similar bound, Qem ~ ^-^ ^ 10"^'^ h{3p — 1) eV cm^'^ sec~^, 
which leaves the bound on rj and mx practically unchanged. Furthermore, constraints from limits on 
CMB distortions and light element abundances from '^He-photodisintegration are comparable to the 
bound from the directly observed diffuse GeV 7-rays [^. That these crude normalizations lead to 
values of r] in the right range suggests that defect models require less fine tuning than decay rates in 
scenarios of metastable massive dark matter. 

As discussed above, in TD scenarios most of the energy is released in the form of EM particles 
and neutrinos. If the X particles decay into a quark and a lepton, the quark hadronizes mostly into 
pions and the ratio of energy release into the neutrino versus EM channel is r ~ 0.3. The energy 
fluence in neutrinos and 7— rays is thus comparable. However, whereas the photons are recycled down 
to the GeV range where their flux is constrained by the EGRET measurement, the neutrino flux is 
practically not changed during propagation and thus reflects the injection spectrum. Its predicted 
level is consistent with all existing upper limits (compare Fig. |2| with Fig. 0) but should be detectable 
by several experiments under construction or in the proposal stage (see Fig. ^). This would allow to 
directly see the quark fragmentation spectrum. 



5 Conclusions 

Ultra-high energy cosmic rays have the potential to open a window to and act as probes of new 
particle physics beyond the Standard Model as well as processes occuring in the early Universe at 
energies close to the Grand Unification scale. Even if their origin will turn out to be attributable to 
astrophysical shock acceleration with no new physics involved, they will still be witnesses of one of 
the most energetic processes in the Universe. The future appears promising and exciting due to the 
anticipated arrival of several large scale experiments. 
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Figure 1: (from work in Ref. |^0[) Cosmogenic neutrino flux per flavor (thick line, assuming maximal 



mixing among all flavors) from primary proton flux (thin line) fitted to the AGASA cosmic ray data |0 
above 3 x 10^^ eV (error bars). The UHECR sources were assumed to inject a E""^ proton spectrum up 
to 10^^ eV with luminosity oc (1 + z)^ up to 2; = 2. Also shown are existing upper limits on the diffuse 
neutrino fluxes from AMANDA AGASA the Fly's Eye [|9] and RICE |0[ experiments, and the 
limit obtained with the Goldstone radio telescope (GLUE) 
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Figure 2: (from work in Ref. [^) Predictions for the differential fluxes of 7— rays (dotted line), nucleons 
(thin solid line), and neutrinos per flavor (thick solid line, assuming maximal mixing among all flavors) 
in a TD model characterized hj p = 1, mx = 2 x 10^^ GeV, and the decay mode X ^ q + q, assuming 
the QCD fragmentation function in MLLA approximation 0], with a fraction of 10% nucleons. The 
calculation used the code described in Ref. |T8| and assumed the minimal URB version consistent with 
observations [26| and an EGMF of 10~^^ G. Cosmic ray data are as in Fig. |l| and the EGRET data on the 
left margin represents the diffuse 7— ray flux between 30 MeV and 100 GeV ||2^. Also shown are expected 
sensitivities of the Auger project currently in construction to electron/muon and tau-neutrinos |5^, and 
the planned projects telescope array |5^, the water-based ANTARES [|^], the ice-based ICECUBE [^9] , 

(we take the latter as 



the fluorescence/ Cerenkov detector MOUNT 
representative also for EUSO 



60[| , and the space based OWL 
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